Backward diode composed of a metallic and semiconducting nanotube 
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The conditions necessary for a nanotube junction connecting a metallic and semiconducting 
nanotube to rectify the current are theoretically investigated. A tight binding model is used for the 
analysis, which includes the Hartree-Fock approximation and the Green's function method. It is 
found that the junction has a behavior similar to the backward diode if the gate electrode is located 
nearby and the Fermi level of the semiconducting tube is near the gap. Such a junction would be 
advantageous since the required length for the rectification could be reduced. 
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Carbon nanotubes (NTsp are promising components 
for nano-electronics because-ithey can either be semicon- 
ducting (S) or metallic (M)0. This has been explaced in 
various experiments, for example, NT transistorsti, NT 
diodeffl, and NT junctionsErQ. The NT junctions made 
by a pentagon-heptagon defect pair can be classified into 
three types: MM, MS, and SS: M and S stand for a metal- 
lic NT and a semiconducting NT, respectively. It is re- 
markable that all the chemical bonds of NT junctionS|^Ee 
essentially sp^ orbitals of the carbon with no impurityErEI, 
which favors a coherent electron flow through the junc- 
tion. For MM-NT junctions, a scaling law of the elec- 
tronic transport was foundcl. On the other hand, MS- 
NT junctions are expected to cause the Schottky barrier 
which is important to rectify the currentlij. In this work, 
we focus on MS-NT junctions connecting the (9,0) tube 
and the (8,0) tube shown in Figj^. Though iii|real exper- 
iments, these junctions are sometimes bentcl, we inves- 
tigate straight junctions because the difference does ipjb 
alter the results qualitatively as was explained in Ref.t3. 

For NT devices, the gate electrode is quite important 
for controlling the transport properties, as was shown 
in an experinpnt where an AFM tip was used as the 
gate electrodelHl. It was recently pointed out that the 
length of the NT junction needs to be larger than the 
nanoscale in order to rectify the current becauae-of the 
large depletion width in the underdoped casdljila. An 
important prospect is to understand whether or not the 
depletion width can be made narrower by approaching 
the gate closer to the NT junction. In this paper, we 
demonstrate that this could indeed be possible because 
the gate electrode screens the Coulomb interaction in the 
junction. The main purpose of this paper is to support 
the validity of this method. 

The current flowing through MS-NT junctions can 
be quite sensitive to the details of the Schottky bar- 
rier caused by the Coulomb interaction. Though 
MS-NT junctians have been recently investigated 
theoreticallytjili the Coulomb interaction was not con- 
sidepid in Ref.113 and the semi-classical model used in 
Ref.tj is not appropriate when the spatial range of 
the potential variation is short. For these reasons we 
have adopted a tight binding model with tt orbitals and 
the Hartree-Fock approximationll3. The corresponding 
Hamiltonian H is represented with the density matrix 
P<j by 

Hi,] = Uj - ^U(\fi -~ rj\)pij + ^ U{\n - rk\){Pk,k - 

k 

(1) 

where fj is the position of the carbon atom j which can 
be determined approximately by the condition that all 
bonds must essentially have the same length. This is jus- 
tified since the detailed atomic structure do not modify 
significantly the electronic structurdj. The first term in 
Eq. (|l|) has a nonzero constant value —t<0 when i and 



j are nearest neighbors and is zero otherwise. The sec- 
ond and third term are the exchange and Hartree term, 
respectively. The Coulomb interaction U(r) is assumed 
to be U{r) = t/ y/l + 4(r/a)^ where a is the bond length. 

The Green's function-method is used for calculating 
the transmission rateESE3. The system is divided into 
three regions: the region near the defects (C), the (8,0) 
tube far from the defects (S), and the (9,0) tube far from 
the defects (M). The corresponding Hamiltonian H is 
written as H = He + Hs + Hm + ts + tM, where tp 
{ p = S 01 M) represents the off-diagonal elements be- 
tween the C region and the p region. For each region 
p, the matrix element (Hp + tp)i,j is written as Tp when 
i and j are nearest neighbors, ep when i = j, and zero 
otherwise. For each iteration step, the parameter ep and 
Tp are determined by the condition that the near and far 
region are connected smoothly, i.e. and Tp are taken 
to be the same as the corresponding averaged value of 
{Hc)i.j with the atom j directly connected to the region 
p. The retarded Green's function G{E) = (E + iS — H)^^ 
can then be represented by 

Gk,i = {E - He - As - AM)k} , (2) 

where k and / belong to the region C. The effect of 
the region S and M are included in the 'self energies', 
Ap = ip{E + i5 - Hp)^^ip. From Eq.(^, the density 
matrix p in the region C can be calculated as 

P=- GT2G*dE+- / G{T2+Ti)G*dE . 

(3) 

Here Tp = z(Ap-A*), Epi = unn{EFMi Eps) and Ef2 = 
mayilEpM , Eps)- Thus, H and p can be calculated self- 
consistent ly with Eqs.(0), (||), and (^). It is worth noting 
that the integrand of the second term in Eq,^ equals the 
DOS because G{Ts + Tm)G* = i(G - G* fK 

Although the electron density pi^i in Eq.(||) is only de- 
fined for the region C, the Hartree terms caused by the 
charges outside this region should be included in He- 
This can be done by repeating periodically that NT unit 
cell in the region C which is nearest to the far region. Fur- 
thermore the electrostatic potential caused by the gate 
electrode can be obtained by the image charge method. 
Here we can choose the electrostatic potential to be zero 
at the gate surface without loss of generality because the 
^^difference, not the absolute value, of the electrostatic po- 
tential is relevant quantity. Note that the nanotubes are 
not neutral by themselves, but the total system includ- 
ing the image charge is neutral. In our model, the gate 
is planar and parallel to the tube axis. Its distance from 
the tube axis is chosen to be three times the radius of 
the (9,0) tubCp^This is much smaller than the gate sep- 
aration in RefJlJ. At this distance, we have verified that 
the orientation of the defects does not change our main 
results. In the following, the defects are assumed to face 
the gate (see Fig.Q). 
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When self-consistency has been achieved, the current 
I is calculated by 

/ = (2e//i) i:T[TsGVMG*]dE . (4) 

The integrand of Eq. (Q) equals the total transmission rate 
T{E) = - Ti,j (E) where Tij is the transmission rate 
from the ith channel to the jth channel. In the energy 
range considered here, both tubes have two channels ex- 
cept for the gap region of the (8,0) tube. Note that T{E) 
may be larger than one, though Tij < 1. 

The electron potential at each atom site is shown by 
circles in Fig.||. The Fermi levels (EpM^Eps) shown by 
straight 

broken lines are (-0.230t, -0.344i), (-0.480t, -0.344i), 
(-0.370i, -0.4800 and (-0.600i, -0.480t) in (a), (b), (c) 
and (d), respectively. These four situations are referred 
to(a),(b),(c) and (d), hereafter. The horizontal axis is 
the position along the tube axis and its range is as same 
as the length of the region C. In Fig.||, the local gap 
energy region (hatched area) and the valence band edge 
(thin curved broken line) in the (8,0) tube are determined 
by the averaged electron potential and the gap width of 
Hs- The Fermi levels correspond to hole doping which 
is not due to dopant atoms, but caused by the source, 
drain and gate electrodes. The regions occupied by the 
electrons are shown by gray area. Though the energy 
range between Eps and EpM is partially occupied (this 
corresponds to the first term of Eq.(|^)), it is not shown 
in Fig.^ explicitly. One may notice that several circles 
show a sharp increase or rapid drop near the interface 
in Fig.^j. These abrupt changes are caused by an excess 
of the electtpn on the pentagon and a deficiency on the 
heptagontlO. 

The origin of the band bending seen in Fig|| (b) and 
(d) is explained qualitatively with the 'shifted Fermi 
level', Epp = Epp — Cp, as follows. Using the en- 
ergy integral of the DOS, the hole density in the far 
regions can be approximated by p'^j ^ Cm\Epm\ and 

~ Cs\l E%S - A2 with the half gap width A. Here 
Cm — Cs because the two tubes have similar radii. Then 
approximate values of the ratio p%/ p\i are 0.32 and 0.44 
for Fig.||(b) and (d), respectively, so p^ < p\j in the 
both situations. The charge distribution connecting p^j 
to p'g causes the electron potential to decrease toward 
the metallic tube, and is responsible for the band bend- 
ing. Comparing Fig.^b) and (d), one can see that both 
the height and the width of the barrier decrease as the 
hole density increases. Qualitatively, this relation is 
similar to that found in Ref.t2l. 

The inset of Fig.^ shows the dependence of the current 
/ on the bias voltage Vsd = {Epu — Eps)/\e\. The six 
open triangles and the five closed circles represent the 
data when Eps = -0.344^ and Eps = -0.480i, respec- 
tively. The data corresponding to Fig.^ are indicated by 
the arrows. The horizontal separation of the triangles or 



the circles is not regular because each Fermi level Epp 
was calculated with a fixed shifted Fermi level Epp. This 
was done since the constant Epp keeps the Pp unchanged 
in the iteration so that it is advantageous to achieve self- 
consistency. 

The open triangles in the inset of Fig.|| indicate that 
the current with constant Eps — — 0.344t is rectified, i.e. 
the resistance is much higher when Vad < than when 
Vsd > 0. For example, Vsd /I ~ 32 kfi for the case (b) 
and Vsd/ 1 ~ 1300 kf7 for the case (a). A positive Vsd: 
however, corresponds to a reverse bias for the thermo- 
electron current in conventional MS junctions, i.e. the 
barrier becomes higher as Vsd increases as can be seen 
in Fig.^. The mechanism of this rewersetL rectification 
is similar to that of the backward diodetS, which can 
be understood by examining the total transmission rate 
T{E) as a function of the energy (Fig.||). Since the nec- 
essary energy range to calculate / with Eq. (^ is between 
EpM and Eps, lines in Fig.^ are limited to this range. 
T{E) in the case (b) is considerably larger for the whole 
bias window because the spatial thickness of the barrier 
is reduced due to the proximity of the gate (see Fig.|^ 
(b)). This leads to a large current for positive Vsd- If the 
barrier had been much thicker as in Ref.E2l, the current 
would have been close to zero. On the other hand, T{E) 
in the case (a) becomes zero for most of the bias window 
since the energy gap of the (8,0) tube blocks the current 
as can be seen in Fig|2| (a) . This example makes it clear 
that the necessary condition for keeping the current small 
when Vsd is negative is to fix Eps near the energy gap. 
Since Eps — — 0.344f is closer to the energy gap than 
Eps = — 0.480i, the triangles show much smaller current 
than the circles in the inset when Vsd < 0. This means 
that a better rectification has been achieved in the for- 
mer case. Here we should mention that —{Eps+Epm)/'2 
waSjfijxed , i.e. Eps depends on Vsd, in each I-V curve in 
Rcf.tj. This is not suitable for the rectification presented 
here. 

In the present paper, we have to use a relatively smaller 
Coulomb interaction ?7(ri--j= t/ \/i + iir/a)"^ than that 
estimated in other paperso, since the long spatial range 
of U (r) makes the convergence difficult. U (r) is probably 
larger in the actual experiments than that used here, but 
it is expected that values of larger U{r) favor our conclu- 
sion. In fact, it is well known that the screening length W 
is determined by the relation W x 1/ y^U{0)Np where 
Np is the DOS at the Fermi level. Though this relation 
itself might not hold for MS-NT junctions, a decrease of 
W with a simultaneous increase of U{0) will hold and 
favor our conclusion. 

In summary, we have calculated the I-V characteristics 
of a typical MS-NT junction. It is found that approach- 
ing the gate electrode to the junction and fixing Eps near 
the energy gap could be an effective metliad-to overcome 
the problem of the large depletion lengthEj'113. The recti- 
fying MS-NT junction discussed here is analogous to the 
backward diode. Though density functional calculations 
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will ultimately be necessary, the results obtained here are 
worth testing in experiments. 
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FIG. 1. The junction connecting the (9,0) tube and the 
(8,0) tube. Above right: setup of the electronic circuit. 

FIG. 2. The electron potential energy at each atom 
site. The Formi levels {Efm,Efs) shown by straight 
broken lines are (-0.230t, -0.344t), (-0.480t, -0.344t), 
(-0.370t, -0.480t) and (-0.600t, -0.480t) in (a), (b), (c) and 
(d), respectively. The unit of the vertical axis is the transfer 
integral t. The horizontal axis represents the position along 
the tube axis and its unit is the lattice constant of graphite, 
i.e. \/3 times the bond length a. Its origin is the interface 
between the (9,0) tube and the (8,0) tube and its range is 
the same as that of the neax region C. The horizontal arrows 
indicate the direction of the electron's flow. The local gap en- 
ergy region, the valence band edge, and the region occupied 
by electrons are shown by the thin curved broken line, the 
hatched area and the gray area, respectively. 

FIG. 3. Total transmission rate as a function of energy 
in the bicis window. The unit of the horizontal cixis is the 
transfer integral t. Inset: I-V characteristics of the junction 
when Efs is fixed to — 0.344f (open triangles) and — 0.480t 
(closed circles). The units of the iiorizoiital and vertical axis 
are 2et/h ~ 0.21 mA and t/\e\ ~ 2.7 V, where e and h arc 
the electron charge and Planck's constant, respectively. Data 
corresponding to Fig.2 are indicated by the arrows. 
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